Tyurin parameters and elliptic analogue of 
nonlinear Schrodinger hierarchy 



Kanehisa Takasaki 
Graduate School of Human and Environmental Studies, Kyoto University 
Yoshida, Sakyo, Kyoto 606-8501, Japan 
E-mail: takasaki@math.li.kyoto-u.ac.jp 

Abstract 

Two "elliptic analogues" of the nonlinear Schrodinger hiererchy are constructed, 
and their status in the Grassmannian perspective of soliton equations is elucidated. 
In addition to the usual fields u, v, these elliptic analogues have new dynamical 
variables called "Tyurin parameters," which are connected with a family of vector 
bundles over the elliptic curve in consideration. The zero-curvature equations of 
these systems are formulated by a sequence of 2 x 2 matrices A n (z), n = 1, 2, . . ., of 
elliptic functions. In addition to a fixed pole at z = 0, these matrices have several 
extra poles. Tyurin parameters consist of the coordinates of those poles and some 
additional parameters that describe the structure of A n (z) , s. Two distinct solutions 
of the auxiliary linear equations are constructed, and shown to form a Riemann- 
Hilbert pair with degeneration points. The Riemann-Hilbert pair is used to define a 
mapping to an infinite dimensional Grassmann variety. The elliptic analogues of the 
nonlinear Schrodinger hierarchy are thereby mapped to a simple dynamical system 
on a special subset of the Grassmann variety. 
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1 Introduction 



Many integrable systems are expressed in the form of a Lax equation dtA(X) = [B(X), A(X)] 
or a zero-curvature equation [d x — A(X), dt — B(X)} = 0, where A(X) and B(X) are ma- 
trices of rational functions of the spectral parameter A. In other words, these Lax or 
zero-curvature equations are defined on the Riemann sphere. Some integrable systems, 
such as the elliptic Calogero-Moser system and the Landau-Lifshitz equation, have a Lax 
or zero-curvature representation defined on a torus, i.e., a complex elliptic curve. One 
will naturally expect to find a generalization to a curve of higher genus. Unfortunately, it 
is well known that such a naive attempt will be confronted with a serious difficulty that 
stems from the Riemann- Roch theorem [TTj . 

Recently, Krichever presented a general scheme for constructing a Lax or zero- curvature 
equation on an algebraic curve T of arbitrary genus $]. A central idea is to allow the 
matrices A, B to have extra "movable" poles 71, ... , 7 r9 G T, where r is the size of the ma- 
trices. Moreover, the matrices A, B at these poles are assumed to have a special structure. 
A set of additional parameters are introduced to parametrize this special structure. The 
coordinates of poles and these parameters are called "Tyurin parameters." This notion 
originates in algebraic geometry of holomorphic vector bundles over algebraic curves JS] , 
and was applied by Krichever and Novikov in 1970's to the study of commutative rings 
of differential operators [3J U\ ■ The aforementioned difficulty can be resolved by adding 
Tyurin parameters as new dynamical variables. 

Once applied to zero-curvature equations, Krichever's method yileds a large class of 
1 + 1 dimensional integrable PDE's. These equations are to be called "soliton equations" 
associated with an algebraic curve (though it is not known whether these equations do have 
a soliton or soliton-like solution). For instance, Krichever illustrates his construction for 
the case of a "field analogue" of the elliptic Calogero-Moser system. This raises a natural 
question: What is the status of these new equations in the Grassmannian perspective of 
soliton equations due to Sato [H] and Segal and Wilson [To] ? 

We address this problem in a simplified setting, namely, zero-curvature equations 
of 2 x 2 matrices defined on an elliptic curve. This system is an analogue of the usual 
nonlinear Schrodinger hierarchy. More precisely, we construct two distinct versions of this 
"elliptic analogue," one being based on Krichever's idea, and the other inspired by the 
work of Enriquez and Rubtsov 2\. Whereas Krichever's construction requires all Tyurin 
parameters to be dynamical variables, Enriquez and Rubtsov keeps the position of poles 
constant and use the other parameters as dynamical variables. In this respect, the elliptic 
analogue a la Enriquez and Rubtsov's is much closer to usual soliton equations. 
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Our strategy is, firstly, to derive a kind of Riemann-Hilbert problem for these systems, 
and secondly, to translate it to the language of an infinite dimensional Grassmann vari- 
ety. This is indeed the procedure that has been used in the literature for many soliton 
equations and some higher dimensional systems; see, e.g., the book of Mason and Wood- 
house [O]. The usual Riemann-Hilbert problem, however, does not work literally in the 
present situation. Whereas the usual Riemann-Hilbert problem is based on triviality of a 
holomorphic vector bundle over the Riemann sphere, the systems formulated by Tyurin 
parameters are obviously related to a nontrivial holomorphic vector bundle over an al- 
gebraic curve of positive genus. An answer to this puzzle can be found in the work of 
Krichever and Novikov E] cited above. They consider a Riemann-Hilbert problem 
with degeneration points] Tyurin parameters are nothing but the geometric data of those 
points. The next task is, therefore, to connect this kind of Riemann-Hilbert problems 
with an infinite dimensional Grassmann variety. Fortunately, a related issue has been 
investigated by Previato and Wilson [13J. They demonstrate therein a Grassmannian 
version of the "dressing method" — a classical technique in soliton theory — to solve a 
Riemann-Hilbert problem of the same type. Moreover, their paper shows what should 
be the "vacuum" (to be "dressed") that corresponds to a holomorphic vector bundle in 
the Tyurin parametrization. Our goal is to develop a similar machinery for the present 
setting. 

This paper is organized as follows. Section 2 is a brief review of the usual nonlinear 
Schrodinger hierarchy. This will serves as a prototype of the subsequent construction. 
Section 3 is devoted to the construction of the first version, a la Krichever, of the elliptic 
analogues. A technical clue is a generating function U(z), which has been used for the 
usual nonlinear Schrodinger hierarchy as well. This enables one to formulate the gener- 
ators of time evolutions systematically. Section 4 deals with an auxiliary linear system 
of the hierarchy and a pair of solutions thereof. This pair of solutions turns out to sat- 
isfy a Riemann-Hilbert problem with degeneration points on the elliptic curve. Section 
5 presents main results of this paper, namely, a Grassmannian perspective of the elliptic 
analogue of the nonlinear Schrodinger hierarchy. An infinite dimensional Grassmann vari- 
ety Gr, a special basepoint ( "vacuum" ) Wo G Gr and the set M. C Gr of "dressed vacua" 
are introduced. The Riemann-Hilbert pair determines a point of Ai, whose motion turns 
out to obey a simple exponential law. The elliptic nonlinear Schrodinger hierarchy is thus 
mapped to a dynamical system on M.. In Section 6, the same story is repeated for the 
elliptic analogue a la Enriquez and Rubtsov. Our conclusion is shown in Section 7. 
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2 Nonlinear Schrodinger hierarchy 



As a prototype of the elliptic analogue, we here review a standard construction of the 
nonlinear Schrodinger hierarchy. Generalities and backgrounds of this kind of construction 
of soliton equations can be found in Frenkel's lectures [3]. 

2.1 ^-matrix 

The construction starts from the A-matrix 

with a rational spectral parameter A G P 1 ; u and v are fields on the x space. In view 
of the homogeneous grading of an underlying loop algebra, it is natural to express this 
matrix as 

A(X) = JX + A {1) (2.2) 

where 

\0 -1 J \ v 

2.2 Generating functions 

A clue of the construction of the hierarchy is a Laurent series 

oo 

U{\) = Un\~ n , U = J, 
n=0 

that satisfies the differential equation 

[d x -A{\), U(X)}=0. (2.3) 

Although this equation itself does not determine U(X) uniquely, there is a good or canonical 
solution that takes the form 

U(X) = 0(A) J0(A)-\ (2.4) 
where 0(A) is a Laurent series of the form 

oo 

0(A) = /+$> n A-" 

n=l 
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and satisfies the differential equations 



40(A) =A(A)0(A)- 0(A) JA. (2.5) 

A solution of (J2.3j) of this form is indeed "good" or "canonical" in the sense that the 
coefficients U n can be calculated from A(X) by a purely algebraic procedure (namely, 
without actually solving differential equations) as follows. Expanded in powers of A, (J2.3j) 
becomes a system of differential equations 

d x U n = JU n+1 - U n+1 J + [A®, U n ] 

for the coefficients U n . On the other hand, if U(X) is written as 1)2.4)1 . the algebraic 
constraint 

U(X) 2 = I (2.6) 
is automatically satisfied. This yields the algebraic relations 

n 

= JU n+1 + U n+1 J + U mUn +l-m 

m=l 

of Un's. One can use these relations to eliminate the term U n+ \J on the right hand sides 
of the foregoing differential equations. The outcome are the recurrence relations 



2JU n+1 = d x U n - [A®,U n ] - £ U m U n+1 - m (2.7) 
that determine U^s successively as 









: 




, u 2 =( 









m=l 



— \uv \u a 



\v x \uv 



etc. Note that the matrix elements of all U n J a thus turn out to be "local" quantities, 
namely, polynomials of ^-derivatives of u and v. 

The coefficients of 0(A) are "nonlocal." To construct 0(A) from A(X), one expands 
()2.5j) to the differential equations 

dAn = [J,0„+i] + A«0„ (2.8) 

for the coefficients and solves them step by step. Actually, this is not so straightforward; 
one has to split 0„ into the diagonal and off-diagonal parts, 

0n = (0n)diag + (0n)ofT-diag) 
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and consider them separately. The differential equation for the coefficients of 0(A) are 
thereby decomposed to the two equations 

^(0„)diag = (^ (1 Vn)diag (2.9) 

and 

<9x(0ji-l)ofr-diag — [J. 

(0n)ofT— diag 

(For convenience, the index n in the second equation has been shifted.) The first equation 
determines (0 n )diag, up to integration constants, if 0i, . . . , 4> n -\ and (0n)ofr-dia g are given. 
The second equation is rather an algebraic equation that determines (0n) o ff-diag from 
0i, • • • , 4> n -i. To construct a solution, therefore, one has to use these equations in a cyclic 
way: 

1. Solve (12.101) for (0 n ) o fr-dia g - 

2. Solve O for (0 n ) diag . 

3. Increase n by 1 and return to step 1. 

The first step of this cycle is to construct (0i) o fr-diag as a solution of (j2.1()|) (n = 1); note 
that the only data necessary here is 0o = /• Starting with this step, one can proceed as 
(0i)off-dia g -> (0i)dia g -> (02)ofT-dia g -> (02)dia g Changing integration constants in 

the solution of ()2.9j) amounts to the right action 0(A) — > 0(A)C(A) by a diagonal matrix 
C(A) = diag(ci(A), 02(A)) of Laurent series with constant coefficients. 

2.3 Construction of hierarchy 

Having constructed the generating function U(X), one can formulate the hierarchy as the 
system of the Lax equations 

[d tn -A n (X), 17(A)] =0, (2.11) 

where A n (X) denotes the "polynomial part" of U(X)X n : 

A n (X) = U Q X n + U^' 1 + ... + U n . (2.12) 

Since U\ = A^ 1 ', Ai(X) coincides with A(X), so that x can be identified with the first time 
variable t\. As we shall show later in a more complicated situation, one can derive the 
zero-curvature equations 

[d tm -A m (X), d tn -A n (X)} = (2.13) 
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from these Lax equations of U(X). Actually, another set of zero-curvature equations, i.e., 

[dt m -A-(X), d tn -A-(\)]=0, (2.14) 

can be derived for the Laurent "tail" 

A-(X) = A n (X) - U(X)X n = -U^X' 1 - U n+2 X~ 2 (2.15) 

as well. These "dual" zero-curvature equations are the Frobenius integrability condition 
of the linear system 

^0(A) = ^(A)0(A), (2.16) 

which thereby determine the time evolutions of 0(A). This linear system turns out to be 
equivalent to the usual auxiliary linear system 

d tn ^(X) = A n (X)^(X) (2.17) 

upon identifying 

oo 

^(A) = 0(A)exp(]Tt„JA") (h = x). 

71=1 

3 Construction of elliptic analogue a la Krichever 

Let T be a nonsingular elliptic curve realized as the torus C/(2a>iZ + 2cj 3 Z), and z the 
complex coordinate of C, which is also understood as a local coordinate of T. The 
polynomial matrices A(X),A n (X) in the nonlinear Schrodinger hierarchy are replaced by 
matrices A(z),A n (z) of meromorphic functions on T. They have a fixed pole at z = 
(which amounts to A = oo in the nonlinear Schrodinger hierarchy) and two "movable" 
poles at z = 71, 72, 71 7^ 72- 

3.1 A-Matrix on elliptic curve 

The role of the A-matrix in the usual nonlinear Schrodinger hierarchy is now played by a 
2x2 matrix A(z) (z G T) of meromorphic functions on T with the following properties: 

1. A(z) has poles at z = 0,71,72 and is holomorphic at other points. 

2. As z^ 0, 

A(z)=( Z ~ 1 "A+Oiz). (3.1) 
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3. As z — > 7 S , s = 1, 2, 



A(z) = ^-^i + 0(l), (3.2) 
^ - 7s 



where o: s and f3 s are two-dimensional column vectors that do not depend on z. ct s 
is normalized as a s = t (a s , 1). 

7 S and a s in this definition are the Tyurin parameters in the present setting, u and v 
are counterparts of those in the nonlinear Schrodinger hierarchy. All these parameters 
are understood to be dynamical, i.e., a function of x (and the time variables t n to be 
introduced later). We have thus altogether six dynamical variables 71, 72, cci, ct2, u, v. 

Lemma 1 If ol\ 7^ «2 ; a matrix A[z) of meromorphic functions onT with these properties 
does exists. It is unique and can be written explicitly in terms of the Weierstrass zeta 
function ((z) as 

A{z) = £ /3s * a .(C(z - 7s) + C(7s)) + ( C(Z) "A (3.3) 

s=l,2 V V ~^ Z > J 

where 

ft = -!-f - 1 V A = -!-f 1 |. (3.4) 

CKi — a 2 V — a 2 / «i - «2 \ «l 

Proof. The defining properties of A(z) imply that -A(z) can be written as 

A(z)= £ /3 S *a.C(* - 7.) + ^) + C, 

s=l,2 

where C is a constant matrix. By the residue theorem, the coefficients have to satisfy the 
linear relation 

£ /3/a a + J = 

s=l,2 

that ensures that -A(z) is single valued on T. Solving these equations for (3 S leads to the 
formula stated in the lemma. On the other hand, matching with the Laurent expansion 
of A(z) at z = yields to the relation 

A« =J2Ps *«sC(-7s) + c, 

s=l,2 

which determines C. 
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The Tyurin parameters 7 S and a s are required to satisfy the equations 

c^ + Tr/3, = (3.5) 
d x l a s + t a s A^ = n s *a„ (3.6) 



where A^' 1 ) stands for the constant term of the Laurent expansion of A(z) at z — j s , 

A^ = lim ( A{z) - S^t) , 
and k s is a constant to be determined by the equation itself. More explicitly, we have 

»x7i = , ^72 = , (3.7) 

d x ai = — 2a 1 ( 12 — d + a\u, d x a 2 = 2a 2 (i2 — v + a%u, (3.8) 

where 

C12 = C(7i) - C(72) - C(7i - 72), 
and the constants n s take the form 

k s = C12 + a s u. (3.9) 

Oil — Oi2 

As Krichever's lemma jHl Lemma 5.2] shows, these equations ensure that the auxiliary 
linear system d x ip(z) = A(z)ip(z) has a 2 x 2 matrix solution that is holomorphic at z = 7 S 
and invertible except at these points. One will notice from ()3.7|) and ()3.8|) that not all of 
the six dynamical variables 71, 72, ai, a 2 , u, v are independent; for instance, one can solve 
()3.8|) for u and v to eliminate u and v as auxiliary dynamical variables. In the following, 
however, we shall treat these six variables on a equal footing. 

3.2 Generating functions 

We now proceed to the construction of two generating functions 

00 00 

(f>{ Z ) = I + ]T <\> n z\ U{z) = J + E UnZ*. 
n=l n=l 

The first generating function (/>(z) is a Laurent series that satisfies the differential 
equation 

d x <j>{z) = A(z)(j)(z) - (f)(z) Jz- 1 . (3.10) 
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u 
v 



Here A(z) is understood to be its Laurent expansion 

oo 

A(z) = Jz- 1 + A {n) z n ~ l (3.11) 

n=l 

at z = 0; the first few coefficients of this expansion read 

V ', 

A ( 2 ) = 1 I aip(7i) - a 2 p(72) p(7i)-p(72) 

ai - «2 y a 1 a 2 (p(7i) ~ K72)) «2p(7i) - "1^(72) 

etc. 

Lemma 2 Laurent series solution <f)(z) of A3.1U\) does exist. 

Proof. Expanded in powers of z, ()3.10j) yields the differential equations 

n+l 
m=l 

for the coefficients <p n . One can decompose these equations into the diagonal and off- 
diagonal parts. The diagonal part becomes the equation 

n+l 

^2: (0n)diag ^ ] i^A ^0n+l— m)diag; 
m=l 

which determines the diagonal part (0 n )dia g of n up to integration constants. The off- 
diagonal part gives the algebraic relation 

n+l 



3r(0n)off-diag = [J, (0n+l)off-diag] + ^ (A (m) (f) 

m=l 

The off-diagonal part (0 n +i)ofr-diag of 4> n +i is thus determined from 0i, • • • , <f> n . Q 
The second generating function £/ (2) can be obtained from <p(z) as 

U(z) = <f>(z) J(f>(z)-\ (3.12) 
which satisfies the differential equation 

[d x -A(z), U(z)}=0 } (3.13) 
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and the algebraic constraint 



U{z) 2 = I. (3.14) 

As we have seen in the case of the nonlinear Schrodinger hierarchy, this algebraic constraint 
singles out a unique Laurent series solution of (|3.13|) . and the Laurent coefficients can be 
calculated by a set of recurrence relations. 

Lemma 3 The coefficients U n of U(z) satisfy the recurrence relations 

n+l n 

2JU n+ i = d x U n — ^2 [A^ m \ U n+ i- m ] — ^2 U m U n+ i^ m . (3.15) 

m=l m=l 

Proof. (|3.13jl yields the differential equations 

n+l 

d x U n = JU n+1 - U n+l J + £ \A- (m \ t/n+l-m] 

m=l 

for the coefficients U n . The algebraic constraint U(z) 2 = 1 gives the algebraic relations 

n 

= JU n+ % + U J n +i + 2! U m U n+ x- m . 

m=l 

Combining them, one obtains the recurrence relation. | | 

One can thus calculate U n 's successively from the Laurent coefficients of A(z) as 

Ut = 



u 
v 



\uv \u x \ i I / p(7i) - p(7 2 ) 

1„, lr 



■\v x \uv j at - a 2 \ aia 2 (p(ji) - p(7a)) 



etc. In particular, the matrix elements of all ?7 n 's turn out to be a polynomial of x- 
derivatives of u, v, j s , a s . 

3.3 Construction of hierarchy 

Generators of time evolution are 2x2 matrices A n (z), n = 1,2, of meromorphic 
functions on T with the following properties: 

1. A n (z) has poles at z — 0,71,72 and is holomorphic at other points. 
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2. As z -> 0, 

A n (z) =f/(z)z~ n + 0(z). (3.16) 

3. As z — > 7 S , s = 1, 2, 

A^) = ^^ + 0(1), (3.17) 

z - 7^ 

where /3 is a two-dimensional column vector that does not depend on z. 

Lemma 4 If ol\ ^ a 2 , a matrix A n (z) of meromorphic functions on T with these proper- 
ties does exists. It is unique and can be written explicitly as 

n—l I i \m 

A n (z) = £ Pn,s *a.(C(* - 7.) + C(7.)) + E ^-j-dTCWUn-i-m + U n . (3.18) 

s=l,2 m=0 m - 

The vectors (3 are determined by the linear equation 

Y, f3 n /<x s + 11^ = (3.19) 

5 = 1,2 

that ensures the single-valuedness of A n (z) on T. 

Proof. Repeat the same reasoning as the case of A(z). Q 

Solving the last linear equation, one can eventually find an explicit form of A n (z) . 
For instance, Ai(z) coincides with A(z), and A 2 (z) takes the form 

A 2 (z) = £ P2,s *a,(C(* - Is) + C(7.)) + M*) + Ui((z) + U 2 , (3.20) 

5 = 1,2 

where 

^ 1 = —( Ua2 ), ^ = —{~ ttai )- (3.21) 

OL\ — Oi 2 \ —V I Oil — OL 2 \ V I 

Let us assume the genericity condition 

«i ^ a 2 (3.22) 

throughout the following consideration. We now formulate an elliptic analogue of the 
nonlinear Schrodinger hierarchy as the system of the Lax equations 

[d tn -A n (z), U(z)]=0 (3.23) 
12 



for the generating function U(z) and the differential equations 

d tn 7 s + Tr/3„/a s = 0, (3.24) 
d tn *a„ + t a s A^ r) = K n>s a s (3.25) 

for the Tyurin parameters. Here A^ 8 ' 1 ' denotes the constant term of the Laurent expansion 
of A n (z) at z = j 8 , i.e., 

2->7 s ^ Z-Js J 

and K n) <( is a constant determined by the differential equation itself. As in the case of 
(|3.5I) and ([3.6)1 . the two equations ()3.24|) and ()3.25|) for the Tyurin parameters are the 
necessary and sufficient conditions for the auxiliary linear system d tn ijj(z) = A n (z)ijj(z) to 
have a 2 x 2 matrix solution that is holomorphic at z = 7 S and invertible except at these 
points. 

3.4 Zero-curvature equations 

Commutativity of the time evolutions in £ n 's is by no means obvious from the construction. 
If one can derive the zero-curvature equations for A n (;2i)'s, commutativity of the time evo- 
lutions is an immediate consequence. As it turns out below, however, the zero- curvature 
equations in the present setting possess richer contents. 

Let us first derive a "dual" expression of the curvature components 

F mn {z) = [d tm - A m (z), d tn - A n {z)}. (3.26) 

Let A^(z) denote the the "tail" part in the Laurent expansion of (j3.16|) . Namely, 

A+(z) = A n (z)-U(z)z- n , (3.27) 

which has a Laurent expansion of the form 

At{z) = {At +1) - U n+1 )z + (4T +2 > - U n+2 )z 2 + ■■■, 

where A^ 's denote the coefficients of the Laurent expansion 

oo 

A n {z) = ]T A^z m ~ n 

of A n (z) at z — 0. The Lax equations (j3.23J) of U(z) can be rewritten in the "dual" form 
d tn U(z) = \A + n {z) - U(z)z' n , U(z)] = [At(z), U{z)\. (3.28) 
The curvature components turn out to have a similar dual expression as follows. 
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Lemma 5 If the Lax equations \3. 2 6 J\) are satisfied, the curvature components F mn (z) can 
be written in the dual form 

F mn {z) = [d tm - A + m (z), d tn - A+(z)]. (3.29) 

Proof. Differentiating A m (z) = U(z)z~ m + A^z) by t n and using the Lax equation 
()3.23|) . one has 

d tn A m (z) = [A n (z),U(z)]z- m + d tn AUz) 
= iA+(z),U(z)}z- m + d tn Ai(z), 

and exchanging m and n, 

d tm A n (z) = [A+(z), U(z)]z- n + d tm A+ n {z). 

As for the commutator [A m (z), A n (z)], 

[A m {z),A n {z)\ = [U(z)z- m + A+(z), U(z)z~ n + A + n (z)] 

= [Al(z),U(z)]z- m - [A+(z),U(z)]z~ n + [A+{z),A+(z)\. 

Collecting these pieces yields the dual expression of the curvature. Q 



Lemma 6 If \3.2J± ) and 4,'/. 25]) are satisfied, the zero- curvature equations 

[d^-A^z), d tn -A n (z)]=0 (3.30) 
can be derived from the Lax equations \3. 23\) . 

Proof. The following method of proof originates in the early work of Krichever and 
Novikov [7j. The curvature component F mn (z) is a matrix of meromorphic functions on 
T. Suppose that F mn (z) turns out to satisfy the following conditions: 

1. F mn (z) is holomorphic at all points of V other than possible poles at 71,72- 

2. Asz^ 0, F mn (z) = O(z). 

3. As z — > 7 S , s = 1,2, 

F rm (z) = ^ ta ' +0(1), 
where j3 mn s is a two-dimensional column vector. 
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Such a matrix of function can be expressed as 

F mn {z) = J2 Pmn,s WC(* " 7s) + C(7.))- 
s=l,2 

By the residue theorem, the coefficients satisfy the relation 

ftmn,s a s 0, 



s=l,2 



which, under the the genericity condition (|3.22jl . imply that f3 mn s = 0, hence F mn (z) = 0. 
Thus the proof is reduced to confirming that F mn (z) does have the three properties. The 
first and second properties are now obvious; in particular, the dual expression of F mn (z) 
and the fact that A£(z) = 0(z) imply that F mn (z) = O(z) as z — » 0. What is left is to 
check the third property. To this end, note that 



Fynn (^) 



d tm - - + 0(z - 7s ), d tn - - A ^ + 0(z - 7 , 



as z — > 7s- Expanded to powers of z— 7s , one can readily see, by 1)3.24)) . that the coefficient 
of [z — 7 S ) ~ 2 vanishes. It is also easy to see, by ()3.25|) . that the coefficient of (z — 7 S ) _1 
is a rank-one matrix of the factorized form (3 mn s t ct s . Q 

One can conversely derive the Lax equations ()3.23)) from the zero-curvature equations. 

Lemma 7 The Lax equations \3. 2'J\i can be derived from the zero- curvature equations 

Proof. Substituting A m (z) = A^(z) + U(z)z~ m in the zero-curvature equation yields 

[d tm - A + m {z) - U(z)z- m , d tn - A n (z)} = 0, 

which one can further rewrite as 

[d tn - A n (z), U(z)} = [d tn - A n (z), d tm - Ai(z)}z m . 

Since A^z) = 0(1) and A n (z) = 0(z~ n ) as z — > 0, the right hand side of the last 
equation is 0(z m ~ n ), so that 

[d tn -A n (z), U{z)]=0{z m - n ). 

Letting m — > oo, one obtains the Lax equation ()3.23)l as expected. Q 

We thus eventually arrive at the following conclusion. 
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Theorem 1 As far as \3.24 ) and \3. 2d}) are satisfied, the Lax equations \3. 23\) and the 



zero- curvature equations \3. 3U\) are equivalent. 

Let us conclude the present consideration with a comment on (|3.24|) and (|3.25|) . Al- 
though these equations look somewhat distinct from the other equations, these equations 
themselves are directly related to the zero-curvature equations 

[d tn -A n (z), d x -A(z)] = Q (3.31) 

between A n (z) and A(z). Namely, if ()3.5|) and ()3.6|) are satisfied (this should be understood 
as part of the definition of A(z)), (|3.24j) and ()3.25|) follow from these zero-curvature 
equations. One can indeed derive these equations from the Laurent expansion of the 
left hand side of ()3.31|) at z = j s . In this respect, one may consider the zero-curvature 
equations (|3.31j) as the defining equation of a hierarchy. This is indeed the way Krichever 
formulates a hierarchy. 



4 Riemann-Hilbert problem 

In this and next sections, we encounter various initial value problems with regard to the 
time variables t = (t\,t2, ■ ■ ■), in which t\ is identified with x. For this reason, let us 
make the notations slightly more strict. Namely, we write a t-dependent quantity always 
indicating its t-dependence explicitly as A n (t,z), j s (t), a s (t), etc. Otherwise, a quantity 
is understood to be independent of t. 

4.1 Laurent series solution of auxiliary linear system 

As a consequence of (|3.29j) . we have the "dual" zero-curvature equations 

[d tm -A+(t,z), d tn -A+(t,z)}=0. (4.1) 

These equations are the Frobenius integrability condition of the linear system 

d tn <P(t,z)=A+(t,z)<P(t,z). (4.2) 

One can redefine the generating function <p(t, z) to satisfy these equations as well. 

Theorem 2 Upon being suitably modified, the generating function (f>(t, z) satisfies the 
forgoing linear system or, equivalently, 

d t J(z) = A n (t, z)(f>(t, z) - 0(t, z)Jz~ n . (4.3) 
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In particular, 

oo 

V>(*, z) = <f>{t, z) exp(^ t n Jz~ n ) (h = x) (4.4) 

n=l 

gives a Laurent series solution of the auxiliary linear system 

d t J)(t,z)=A n (t,z)il)(t,z). (4.5) 

Proof. One can construct a Laurent series <p(t, z) — I + + • • • as a solution of the 
initial value problem 

dt n <P{t,z) = A+(t,z)ij>(t,z), (f>(t,z)\ t2=t3 =-=o = (/)(t,z)\ t2=t3= ... =0 . 

The Frobenius integrability condition of this system is ensured by the zero-curvature 
equation of A+(t, z)'s. Moreover, since A£(t, z) = O(z), the solution persists to be of the 
form J + 0(z). Now consider the new Laurent series 

U(t,z)=^t,z)J4>(t,z)-\ 

which satisfies the differential equations 

d tn U(t, z) = [d t J(t, z) ■ z)~\ U{t, z)\ = [A+(t, z), U(t, z)}. 

On the other hand, one knows that U(t, z), too, satisfies differential equations of the same 
form, i.e., (|3.28jl . Since U(t, z) and U(t, z) have the same initial data at t-i — t$ — ■ ■ ■ — 0, 
uniqueness of solution of the initial value problem implies that U(t,z) = U(t,z), i.e., 

U(t, z) = 4>{t,z)J4>{t, z)~\ 

so that one can rewrite the foregoing differential equation for (f)(t, z) as 

d t J(t,z) = (A n (t,z)-U(t,z)z- n )0(t,z) = A n (t,z)0(t,z)-4>(t,z)Jz- n . 

Thus (f)(t, z) turns out to fulfill all requirements. 

4.2 Global solution of auxiliary linear system 

The Laurent series solution ip(t, z) of the auxiliary linear system, by its nature, carries no 
information on the global structure of A n (t, z)'s on T. To fill this gap, we now introduce 
another solution x{t, z) that is globally defined on V with several singular points. As it 
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turns out, these two distinct solutions of the same auxiliary linear system play the role of 
the Riemann-Hilbert (or factorization) pair in the usual nonlinear Schrodinger hierarchy. 

To avoid delicate problems, we assume in the following that the solutions of the hier- 
archy under consideration are (real or complex) analytic in a neighborhood of the initial 
point t — 0. 

x(t,z), by definition, is a solution of the auxiliary linear system 

d tnX (t,z)=A n (t,z) X (t,z) (4.6) 
that satisfies the initial condition 

X (0,z) = I. (4.7) 

Since the auxiliary linear system is a collection of ordinary differential equations, any 
solution remains nonsingular as far as the coefficients of the equations are nonsingular. 
Consequently, if z is in a subset of T where A n (0,z)'s are holomorphic, such a solution 
x(t, z) does exists in a (possibly small) neighborhood of t = in the t-space. Since all 
singularities of A n (0, z) on V are located at the three points 0, 7i(0), 72(0), we can conclude 
that the singularities of x (t, z) on V are confined to a neighborhood of these three points 
as far as t is sufficiently close to 0. 

To elucidate the nature of singularities on T more precisely, we expand x{^i z ) hito a 
Taylor series at t — and examine the Taylor coefficients as a function of z. Note that 
this is reasonable, because this Taylor series has a nonzero radius of convergence as far as 
^0,7i(0), 72 (0). 

The Taylor coefficients of X (t, z) at t — can be evaluated by successively differenti- 
ating the differential equations as 

d t n x(t,z) = A n (t,z)x(t,z), 
dt m d tnX {t,z) = (d tm A n (t,z)+A n (t,z)A m (t,z)) X (t,z), 
dt k d tm d tn x(t, z) = (d tk d tm A n (t,z)+d tk (A n (t,z)A m (t,z)) 

+ (d tm A n (t, z))A k {t, z) + A n (t, z)A m (t, z)A k (t, z)) X (t, z), 

etc. Letting t — 0, we are left with a noncommutative polynomial of derivatives of A n 's. 
We can deduce from these calculations the following precise information. 

Lemma 8 The derivatives dt ni ■ ■ ■ dt np x{t, z)\t=o of all orders ofx{t, z) att = are a ma- 
trix of meromorphic functions of z onT with poles at z = 0, 71 (0), 72(0) and holomorphic 
at other points. As z — > 7 S (0) ; s = 1, 2, 

dt ni ---d tnpX (t,z)\ t = = ^ 1 7^ Q t p (0) +Q(l), (4.8) 
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where (3 n ... n s (0) is a two-dimensional constant column vector. 

Proof. As illustrated above, the derivatives of x(t, z) of all order can be written as 

dt ni ■ ■ ■ d tnpX (t, z) = A nu ... tnp (t, z) X (t, z). (4.9) 
Differentiating this equation by t m yields the recurrence relations 

A m ,ni,---,n p 

(t,z)+A (t, z)A m (t, z). 

for the coefficients A ni ^... tnp (t, z). One can prove, by induction on p, that A nit ...^ np (t : z) is 
a matrix of meromorphic functions of z on T with poles at z = 0, ji(t), 72 (i), and 

^,.,n p (t, «) = ' "'--'r . 1 ^ + 0(1) (4.10) 

Z 7«v7 

as z -> 7s (0 j where /3 nii ... npjS (t) is a two-dimensional column vector. Assume that the 
Laurent expansion ()4.10|) holds for A nit ... inp (t, z). The Laurent expansion of A m ni ... np (t, z) 
can be read off from the recurrence relation as 

A m ,ni,---,n p (t^ Z) 

= /3 nij .., w (*) (d tmls (t) + 'a,®?,®) (z - 7s (t)r 2 

+ (d tm /3 ni; ... inp>s (t) • + /3 ni ,., npiS (t) • d tm WO + 

+ 0(1), 

where A^}, iH (t) denotes the constant term in the Laurent expansion ([4. 10)1 . By (|3.24j) . 
the coefficient of (z — 7 s (t))~ 2 vanishes; by (J3.25|) . the terms containing d tjn t ct s (t) and 
t cx s (t)A^ L ' 1 '(t) in the coefficient of (z — 7s(t)) _1 cancel out. Thus A m>nir __ irip (t, z), too, 
turns out to have a Laurent expansion of the expected form. This completes the proof 
of (|4.10p . Lastly, letting t = in ()4.9|) . one eventually arrives at the statement of the 
lemma, n 

All Taylor coefficients of x(^, z) at t = thus turn out to have poles at the same 
position, namely, the three points 0, 7i(0), 72(0). Moreover, whereas the order of pole at 
z = is unbounded, the poles at z = 7«(0), s = 1,2, are of the first order. Accordingly, 
x{t, z) has an essential singularity at z = and simple poles at the other two points. The 
leading part of the Laurent expansion at z — 7 S (0) takes the familiar form 
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where (3 xs (t) is a two-dimensional column vector that depends on t. Note that the pole 
of x(t, z) at z — 7 S (0) disappears when t — (because %(0, 2) = /). 
Lastly, let us mention another important property of x(t,z). 

Lemma 9 det x(t,z) is a meromorphic function on T with simple poles at z = 7 S (0), 
s = 1,2, and simple zeroes at z = j s (t), s = 1,2. t a s (t) is a left null vector ofx(t, 7s(0)- 

Proof. The auxiliary linear system d tn x(t, z) = A n (t, z)x(t, z) induces the linear system 

d tn det X (t, z) = Tr A n (t, z) det X (t, z) 

for n — 1, 2, Since A n (t, z) = U (t, z)z~ n + 0(z) as z — > and Tr [/(£, = 0, one finds 

that the coefficients of this linear system for det x(t,z) has no singularity at z — 0, but 
rather a zero, namely, 

Tr A n (t,z) = 0(z) (z^O). 

This implies that 2) has no singularity at 2; = 0. In view of the initial condition 
x(0, 2) = /, one can conclude that det z)\ z=0 = 1. One can thus confirm that 
detx(t, 2) is a meromorphic function on T with poles at 2 = 7«(t), s = 1, 2, and holomor- 
phic at other points. Since the residue matrix of x (t, z) at z = 7 s (t) is a rank-one matrix, 
det x(t,z) has a simple pole there. The position of zeroes of det x(t,z) can be deduced 
from the the linear equation 

d x x(t,z) = A{t,z) X (t,z) 

(or from any any member of the auxiliary linear system). Extracting the residue at 
z — Isif) yields the relation 

= {3 s (t) t as(t)x(t, ls (t)) 
which, because (3 s (t) 7^ 0, reduces to the relation 

t a s (t) X (t, ls (t)) = 0. 

Thus t ac s (t) turns out to be a left null vector of x(^7s(0)- O n ^ ne other hand, rewriting 
the linear system as 

A = d xX (t,z) ■ X (t,z)-\ 

one can see that the zeroes 7 s (i) of det X (t, z) are simple. If they are a multiple zero, the 
matrix A will have a multiple pole; this contradicts the construction of the matrix A. [7 
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These results show that x{t, z) is exactly the solution mentioned in Krichever's lemma 
[SJ Lemma 5.2], namely a matrix solution holomorphic at the movable poles of A(t, z). 
In summary, x{t, z) has the following properties. 

Theorem 3 xifi z ) has an essential singularity at z = and simple poles at z = 7 S (0), 
s = 1,2, and is holomorphic at other points of V . As z — > 7 S (0) ; x(^,z) behaves as 
\4-H\ ) shows. Moreover, det x(t, z) is a meromorphic function on T with simple poles at 
z — 7a(0), s — 1, 2, and simple zeros at z = 7 s (t), s — 1, 2. t o: s (t) is a left null vector of 
X(t,ls(t)). 

4.3 Riemann-Hilbert problem with degeneration points 

We now have two distinct solutions of the same linear system, namely, the Laurent series 
solution ip(t,z) and the solution x(t, z) carrying global information on V. The "matrix 
ratio" of these these two solutions is a constant matrix, i.e., 

d tn (x(MrV(M)exp(fXJ*-")) =0. 

Equating this matrix ratio with its value at t = 0, we are led to the relation 

oo 

X (t,z)- 1 <j>{t,z)e 3 v( y £t n Jz- n )=<j>(0,z) (4.12) 

n=l 

or, equivalently, 

oo 

0(0, z) exp(- ]T t n Jz- n ) = X (t, *) _ V(t, z)- (4-13) 

n=l 

The last relation may be thought of as a kind of Riemann-Hilbert problem concerning 
a small circle \z\ = a on the torus T. The input of this problem are the initial values 
7 s (0),a s (0) and 0(0,2;). The left hand side of ()4.13|) is a GL(2, C)-valued function on 
the circle, in other words, a GL(2, C) loop group element. The problem is to factorize it 
to two factors. The second factor 0(t, z) is a loop group element that can be extended 
to a matrix of holomorphic functions on the inside of the circle. The first factor x(t, z) 
is a loop group element that can be similarly extended to the outside of the circle, but 
not holomorphic everywhere; x(t, z ) is required to have poles at z = 7 S (0), s = 1,2, 
with the structure described in (|4.11J) . Moreover, in addition to these poles, x{pi z ) have 
degeneration points, i.e., zeros of the determinant at z — 7 s (t), s — 1, 2. These zeroes are 
nothing but the poles of A n {t, 2;)'s. 
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Thus the Riemann-Hilbert problem relevant to the present setting is a Riemann- 
Hilbert problem with movable degeneration points (and extra fixed poles) on a torus. 
A similar Riemann-Hilbert problem appears in Krichever's work jHj on commutative rings 
of differential operators. In that case, the Riemann-Hilbert problem is formulated on 
the "spectral curve" of the commutative ring under consideration, and the genus of the 
spectral curve can be an arbitrary positive integer. 

Krichever converts the Riemann-Hilbert problem to an integral equation and solves it 
by a standard procedure. The same method can be applied to the present setting as well, 
though we shall not seek this approach here. An alternative approach, as demonstrated 
by Previato and Wilson ^3] , is to translate the Riemann-Hilbert problem to the language 
of an infinite dimensional Grassmann variety. We shall present this method in the next 
section. 

4.4 Back to hierarchy 

It will be instructive to show how to derive a solution of (J3.23)) . ()3.24j) and (jri.25j) from the 
Riemann-Hilbert problem. This is more or less parallel to the procedure that Krichever 
and Novikov employ in their work [SJ[7|. 

Notice, first of all, that x{^i z ) * s a matrix version of the "vector Baker- Akhiezer 
function" in their terminology. This is an immediate consequence of the Riemann-Hilbert 
problem: z) has an essential singularity of the exponential type at z = 0, and fixed 
simple poles at z = 7 S (0), s = 1,2. Accordingly, the determinant det x(t,z) has zeros at 
7 S (£), s = 1,2, that depend on t. Let us consider the generic situation where 7 s (t)'s are 
simple zeros of det x(t, z). The matrices A n (t, z), now defined by 

A n (t, z) = d tn x{t, z) ■ x(t, z)' 1 , 

thereby has simple poles at 7s(£). As simple linear algebraic calculations show, the residue 
of x{t, z)~ l at the degeneration point 7 s (t) is a rank-one matrix. Consequently, the residue 
of A n (t, z), too, is a rank-one matrix and takes the factorized form f3 ns {t) t ct s (t) with a 
common vector a. s (t) independent of n. The dynamical Tyurin parameters 7 s (t), ce s (t), 
s = 1, 2, are thus obtained. According to a general theorem of Krichever and Novikov (re- 
stated in Krichever's recent paper [Hj), these parameters satisfy the differential equations 
(|3~24l and (ET23J) . 

One can now derive the Lax equation ()3.23|) as follows. Differentiating the Riemann- 
Hilbert relation 1)4.13)1 yields another expression of A n (t,z), 

A n (t, z) = d tn (j)(t, z) ■ cj)(t, zY x + U (t, z)z~ n , 
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where U(t, z) is defined as 



U(t,z) = 0(t,z)J0(t,z)- 1 . 

The Lax equations ()3.23|) are thereby satisfied automatically. Moreover, the second ex- 
pression of A n (t, z) also shows the singular behavior of A n (t, z) as z — > 0: 

A n (t,z) =U(t,z)z~ n + 0{z). 

Thus A n (t, z)'s turn out to have all properties that we have assumed in the construction 
of the hierarchy. 

5 Grassmannian perspective 

We here translate the Riemann-Hilbert problem to the language of an infinite dimensional 
Grassmann variety. This leads to a mapping of the elliptic nonlinear Schrodinger hierarchy 
to a multi-time dynamical system on a subset (the set of dressed vacua) of the infinite 
dimensional Grassmann variety. 

5.1 Formulation of Grassmann variety 

Two different models of infinite dimensional Grassmann varieties have been used in the 
literature of integrable systems. One is Sato's algebraic or complex analytic model based 
on a vector space of (formal or convergent) Laurent series The other is Segal and 
Wilson's functional analytic model based on the Hilbert space of square-integrable func- 
tions on a circle [T^j. Which to choose is rather a problem of taste; both of them work 
well in the present context. Let us use Sato's model in the following. Actually, Sato's 
formulation contains a continuous family of different models. Among them, we choose 
one of the presumably simplest models. 

Let V denote the vector space of all 2 x 2 matrices of Laurent series 

oo 

X(z)= ]T X n z n , A n Ggl(2,C), 

n=— oo 

that converges in a neighborhood of z = except at z = 0; gl(2, C) denotes the vector 
space of 2 x 2 complex matrices without any algebraic constraints. This vector space 
is a matrix analogue of V 8 * 111 ^ 00 ' in Sato's list of models [14]; as noted therein, one can 
introduce a natural linear topology in this vector space. 
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We construct an infinite dimensional Grassmann variety Gr from this vector space V 
and the vector subspace 



of all X(z) G V that are holomorphic and vanish at z — 0. The Grassmann variety Gr 
consists of all closed vector subspaces W C V for which the composition of the inclusion 
map W <-^V and the canonical projection V — > V/V + is a Fredholm map of index 0: 



Gr = {W C V | dimKer(iy -> V/V+) = dimCoker(iy -> < oo}. (5.2) 



The so called "big cell" Gr° C Gr is an open subset that consists of subspaces for which 
the map W —> V/V+ is an isomorphism: 



5.2 Vacuum and dressing 

Following the idea of Previato and Wilson 13 J, we now introduce a special element 
Wo(7, at) of the big cell determined by constant Tyurin parameters 7 = (71,72) and 
a = («i, CK2). This is a matrix version of the "vacuum" that Previato and Wilson suggest 
to use for a holomorphic vector bundle in the Tyurin parametrization. 

Lemma 10 Let 7 = (71, 72) be a pair of distinct points ofT, 71 7^ 72, and a = («i, 02) a 
pair of constants satisfying the genericity condition a\ 7^ «2 ■ Then, for any integer n > 
and the matrix indices i,j = 1,2, there is a unique 2x2 matrix w nt ij(z) of meromorphic 
functions on T with the following properties: 

1- w n .ij(z) has poles at z = 0,71,72 and is holomorphic at other points. 

2. As z — > 0, w n> ij(z) = Eijz~ n + O(z), where Eij, i,j = 1, 2, are the standard basis of 



V + = {X(z) e V I X n = for n < 0} 



(5.1) 



Gr° = {W e Gr | W ~ V/V+}. 



(5.3) 



gl(2,C). 



3. As z 



Is, S 



1,2 




n,ij,s s 



+ 0(1) 



z - 7s 



where ct s 



t (a s , I), and (3. 



is another two-dimensional constant column vector. 



The subspace 



W (j,a) = (w niij (z) I n > 0, i,j = 1,2 



(5.4) 



spanned by (the Laurent series of) w n ^j(z) 's is an element of the big cell. 
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Proof. One can confirm the existence and uniqueness of w nt ij(z) in the same way as the 
case of A(z) and A n (z). Since the leading terms EijZ~ n of the Laurent expansion at z = 
are in one-to-one correspondence with the elements of the standard basis {Eijz~ n \ n > 
0, i,j = 1,2} of V/V+, the linear map Wo(y,a) — > V/V+ is obviously surjective. To 
prove the injectivity, note that any element X(z) of Wo (7, a ) H V" + is a matrix of functions 
holomorphic at all points of Y other than possible poles at z — 71,72, behaves as 

x{z) = Px^± + (l) 

at these points (where (5 X s is a two-dimensional column vector), and has a zero at z — 0. 
Such a matrix of function is equal to as one can see by the same reasoning as the proof of 
the zero-curvature equation ()3.30|) . Therefore Wo(7, a) Pi V + — {0}, hence the injectivity 
of the linear map Wo(7, a) — > V/V + follows. rj 

This special base point Wo (7, a) of the big cell plays the role of vacuum in the "dressing 
method." This is a complicated vacuum with nontrivial structure that stems from an 
underlying holomorphic vector bundle over T. We "dress" this vacuum to obtain an 
element W of the big cell that represents a general solution of our hierarchy. Dressing is 
achieved by multiplying a Laurent series (f>(z) from the right side as 

00 

W = W {j,a)<t>{z), (f>(z) = I + £ 4> n z n , n egl(2,C). (5.5) 

n=l 

Our goal in the following is to show that our hierarchy can be mapped to a multi-time 
dynamical system on the set 

M = {W G Gr° I W = W (r/,a)(f)(z), 

7 = (7i, 72) € r 2 , a = (a h a 2 ) £ C 2 , 71 ^ 72 <*i ^ a 2 , n e gl(2, C)} (5.6) 

of these dressed vacua. 

5.3 Interpretation of Riemann-Hilbert problem 

We now translate the Riemann-Hilbert problem (|4.13J) to the language of dressed vacua. 
Because of several reasons, the following consideration is limited to a small neighborhood 
of t = 0. Firstly, this is to ensure that the conditions 71 (t) 7^ 72 (t) and ai(t) 7^ a^t) are 
satisfied; this issue is related to boundaries of the Tyurin parametrization of holomorphic 
vector bundles. Secondly, if t gets large, the dressed vacuum W{t) e M. representing a 
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solution of ()4.13|) can hit the boundary of the big cell, so that more careful analysis is 
required. 

The first step is the following. 

Lemma 11 W (j{t),a(t)) X (t,z) C W o (a(0), 7 (0)). 

Proof. Let w ni ij(t, z), n > 0, i,j = 1, 2, denote the elements of the basis of Wo(j(t), a(t)) 
defined in Lemma ITU1 w n<i j(t,z) has poles at z — 0, 7i(t), 7 2 (t), and behaves as 

W n ,ij{t t Z) = j-r h 0(1) 

as z — > 7a(t). Upon multiplication with 2), the poles at z = J s (t) are cancelled out 
because t ac s (t) is a left null vector of x(t>ls(t)) (see Theorem EJ). Thus one finds that 
Wn,ij(t, z)x(t, z) has an essential singularity at z = 0, simple poles at z = 7 S (0), s = 1,2, 
and is holomorphic at other points of V . The leading part of the Laurent expansion at 
z = 7 S (0) takes the form 

^(MMM)-^ ? - (0)) Y ),O,(0) ^W- 

so that the residue matrix has such a factorized form as (column vector) • t o: s (0). One can 
thus confirm that w nii j(t, z)x{t, z) fulfills all conditions to be an element of ^0(7(0), ct(0)). 

□ 

The next step is to show that the inclusion relation in this lemma is actually an 
equality. To this end, we prove the following lemmas. 

Lemma 12 %(£, z)~ x has an essential singularity at z = 0, simple poles at 7 s (t), s = 1,2, 
and is holomorphic at other points. As z — > 7 S (£), 

y(M) - 1= ^„W'a,W +0(1)| (5?) 
where fl x -i s {t) is a two-dimensional column vector. 

Proof. It is shown in Theorem El that t a s (t) is a left null vector of x(t, 7 s (t)). A clue to 
the proof of the lemma is the fact that the left null space (i.e., the left zero-eigenspace) 
of x(t, 7s(i)) is, actually, one-dimensional and spanned by t oc s (t). If the left null space 
is two-dimensional, x(^7s(0) itself is a zero matrix, so that det z) has a double zero 
at z = 7 S (t); this contradicts the present setting. [Remark: The same reasoning holds for 
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an r x r analogue of the present well. Namely, if the left null space of x{t>ls(t)) 

has k dimensions, then det x{t, z) has a zero of the k-th order at z = 7a(t)-] Bearing this 
fact in mind, one can prove the statement of the lemma as follows. Theorem 0] implies 
that 7a (t) is a simple zero of Extracting the residue from the obvious identity 

z) -1 x(£, z) = I yields the relation 

Res x{t,z)- 1 dz- X {t,ls{t)) = 0. 

z=-y s (t) 

This implies that the residue matrix of x{ki z )~ l at z = 7«(t) is a rank-one matrix of 
the factorized form (column vector) • (row vector). The row vector on the right side is 
accordingly a left null vector of x{t, Js(t))- By the aforementioned fact, one can choose this 
row vector to be equal to t a s (t). Thus the residue matrix turns out to have a factorized 
form as shown in the statement of the lemma. The other properties of x(t, z)~ x , too, can 
be readily derived from Theorem El | | 

Lemma 13 t o: s (0) is a left null vector of x{t, <2) -1 | z=7s (o)- 

Proof. The identity x(t, z )x(t, z)^ 1 = I yields the relation 

Res x(t,z)dz ■ x(t,z)~ l \ z=ls{()) = 0, 

z=7s(0) 

which, by (|4.11|) . takes the form 

^(*) t «*(0)x(M)" 1 U.(o) = 0. 
Since (3 x s (t) 7^ 0, this implies that 

t a a (0)x(*,«)- 1 |*= 7 .(o) = 0. 

□ 

These lemmas show that the inverse matrix x{t, z)~ l has essentially the same prop- 
erties as x(ti z ) except that the position of poles and degeneration points are exchanged. 
Consequently, one can repeat the proof of Lemma[TTJ replacing the role of x(t, z), W (j(t), a(t)) 
and W (7(0),a(0)) with those of x(t,z)- l } W {j{Q), a(0)) and W (r/(t),a(t)), to derive 
the inclusion relation 

W (7(0),a(0))x(t,^)- 1 C Wo{l(t), a(t)). 
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Thus the equality 



W„(7(t), <*{t))x{t, z) = W (a(0), 7 (0)). (5.8) 

follows as expected. 

Having this equality, one can readily convert the Riemann-Hilbert problem to the 
language of dressed vacua as follows. The Riemann-Hilbert relation 1)4.13)1 yields the 
relation 

oo 

W (j(t), a(t))0(t, z) = W (r/(t), a(t)) X (t, z)<f>(0, z) exp(- £ t n Jz' n ). 

n=l 

By (|5.8jl . H / o(7(t), a(t)) absorbs z) to become Wo(7(0), a(0)). The outcome is the 
relation 

oo 

Wo(7(*), «(t))0(t, z) = W (7(0), a(O))0(O, z) exp(- £ t n J*"") , 

n=l 

which means that the dressed vacuum W(t) = Wo(j(t), a(t))(f>(t, z) G M. obeys the 
exponential law 

oo 

W{t) = W(0)exp(-Y / t n Jz- n ). (5.9) 

n=l 

Conversely, one can obtain a solution of the Riemann-Hilbert problem from the ex- 
ponential flows (J5.9)) as follows. (This is a variation of the dressing method of Previato 
and Wilson 13 .) Given a set of initial values 7(0), a(0) and 0(0, z), let us consider the 
exponential flows ()5.9)) sending W(0) = Wo(7(0), a(O))0(O, z) to W(t). A clue is, again, 
the fact that W(t) remains in the big cell as far as t is sufficiently small. In that case, 
the linear map W(t) — > V/V + is an isomorphism. Let <p(t, z) denote the inverse image of 
/ G V/V+ by this isomorphism. Being equal to I modulo V+, (f>(t, z) is a Laurent series of 
the form 

00 

71=1 

On the other hand, as an element of 

00 

W(t) = W o (7(0), a(O))0(O, z) exp(- £ t n J*""), 

n=l 

0(t, 2) can also be written as 

00 

0(t, z) = X (t, z)4>(0, z) exp(- tnJz~ n ) 

n=l 
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with an element x(t, z ) °f Wo(7(0), a(0)). Recalling the definition of Wo (7(0), «(())), one 
finds that z) is a matrix of functions with all properties in the statement of Theorem 
El The associated Tyurin parameters (j s (t), a s (t)) are determined as the position of zeros 
of x{t, z) an d the normalized left null vector of x(A z) at those degeneration points. 

Thus we have been able to show the following fundamental picture of our hierarchy 
as a dynamical system embedded in the Grassmann variety. 

Theorem 4 The elliptic analogue of the nonlinear Schrddinger hierarchy can be mapped, 
by the correspondence W(t) = Woi'jit), a(t))<j)(t, z), to a dynamical system on the set A4 
of dressed vacua in the Grassmann variety Gr. The motion ofW(t) obeys the exponential 
law \5.9\) . Conversely, the exponential flows on M. yield a solution of the Riemann-Hilbert 
problem, \4.13j) . 

Let us conclude this section with a few remarks. 

1. Our approach owes much to the work of Previato and Wilson They use a similar 
Grassmannian version of the dressing method as a tool to reformulate the work 
of Krichever and Novikov El El on commutative rings of differential operators. 
Accordingly, the detail of the dressing procedure is quite different from ours. In 
particular, they take Krichever's "algebraic spectral data" [S] as the input; dressing 
is achieved by a matrix solution of linear differential equations determined by these 
data. In our case, the dressing matrix is the product of 0(0, z) and the exponential 
matrix generating the exponential flows 1)5.9)1 . 

2. In every aspect, the construction of the mapping to the Grassmann variety is related 
to the geometry of holomorphic vector bundles over T. First of all, the Tyurin 
parameters (7 s (t), oe s (t)) themselves correspond to a holomorphic vector bundle that 
deforms as t varies. The subspace Wo (7, a) C V can be identified with the space 
of holomorphic sections of the associated sl(2, C) bundle over the punctured torus 
r\{z = 0}. (j)(t, z) is related to changing local trivialization of this bundle at z — 0. 
Note, in particular, that the data of local trivialization plays the role of dynamical 
variables. This is to be contrasted with the work of Previato and Wilson; in their 
case, a set of functional parameters in the algebraic spectral data play a similar role 
in place of the data of local trivialization. In this respect, our approach is more 
close to Li and Mulase's approach [121 E] to the classification of commutative rings 
of differential operators, in which the choice of local trivialization is treated as an 
independent data. 
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6 Construction a la Enriquez and Rubtsov 



Enriquez and Rubtsov j2] parametrize the sl(2, C) Hitchin system on an algebraic curve of 
genus g > 2 j3j by 3g (rather than 2g) pairs (7 S , a s ), s = 1, . . . , 3g, of Tyurin parameters. 
The roles of parameters are also different from Krichever's formulation. Namely, whereas 
the directional vectors oz s = t (a s , 1) remain dynamical, the poles 7 S are fixed. 

We borrow their idea to construct another elliptic analogue of the nonlinear Schrodinger 
hierarchy. This hierarchy has three pairs (7 s ,a s ), s = 1,2,3, as Tyurin parameters; 7 S 
are constant and a s are variables. In addition to these Tyurin parameters, the hierarchy 
contains the nonlinear Schrodinger fields u, v. The foregoing consideration on the elliptic 
analogue of the Krichever type can be extended to this case with minimal modifications. 



6.1 Construction of ^-matrix 

The v4-matrix A(z) is a 2 x 2 matrix of meromorphic functions on T characterized by the 
following properties: 

1. A(z) has poles at z — 0,71,72,73 and holomorphic at other points. 

2. As z -»• 0, 

A{z)=[ Z ~ X \\+0{z). 
\ v —z I 

3. As z — » 7 S , s — 1, 2, 3, 

.41 -.■) = *. I a \ 1 I c.-'.r' + oai (en 

-ttg — a s 



where X s is a constant to be determined below. 



One can write A(z) itself more explicitly as 

A(z) = E A s ( a \ 1 ) (((z - 7s ) + C(7.)) + [ C{Z) (6-2) 
8=1,2,3 V -°r 8 -a a J \ v -((z) J 

By the residue theorem, the coefficients have to satisfy the linear equations 
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This yields the three linear equations 



^ A s = 0, = _1 > = 0, 

8=1,2,3 s=l,2,3 3=1,2,3 



which can be solve for A.,'s as 
Ax 



2 2 2 2 2 2 

«3 - a 2 x _ «1 - a 3 \ _Ot 2 -a\ 
M — : , ^3 



A "' A 
as far as the Vandermonde determinant 

A = A(ai, a 2 ,a 3 ) = 
does not vanish. This condition 

A ^ 



111 

«1 «2 «3 

2 2 3 
a\ Oi2 «3 



(6.4) 



(6.5) 



(6.6) 



is the "genericity condition" in the present setting. We assume this condition throughout 
the consideration in the following. 

The Tyurin parameters are requirered to satisfy the differential equations (|3.5|) and 
f)3.6)l . Note that the residue matrices of A(z) at z = 7 S can be factorized as 



A 



s I 2 

—a —a s 



X, 



-a. 



a, 1 



so that the role of f3 s is now played by X s *(1, — a s ). (|3.5jl reduces to 

O.,-, = -TrA, ( 1 | =0. (6.7) 

thus being consistent with the assumption that 7 S are understood to be constant. On the 
other hand, ()3.6j) takes the form 

d x a s = Y X r (a s - a r ) 2 (C(7s - 7r) + C(7r)) + «s M ~ 2a sC(ls) + v (6.8) 
with the constant k s uniquely determined as 

K s = M a * ~ ar)(C(ls ~ 7r) + CM) + & S U ~ ({"f s ). (6.9) 
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6.2 Construction of hierarchy 

The construction of time evolutions is fully parallel to the previous case. 
Firstly, we construct a 2 x 2 matrix of generating functions 



U(z) = J2 UnZ n , U = J, 
n=l 

as a solution of the equations 

[d x -A(z), U(z)]=0, U{zf = I. 

The coefficients U n are uniquely determined by a set of recurrence relations; the matrix 
elements thus turn out to be a differential polynomial of a s (s = 1, 2, 3), u and v. 

Having this generating function as local data at z — 0, we now proceed to the construc- 
tion of the generators A n (z) of time evolutions. A n (z) is a 2 x 2 matrix of meromorphic 
functions on T with the following properties: 

1. A n (z) has poles at z = 0,71,72,73 and holomorphic at other points. 

2. Az^O, 

A n (z)=U(z)z- n + 0(z). 

3. As z — > 7 S , s — 1, 2, 3, 

A n (z) = \ n J a \ 1 )( z - 7s )-i + C(l), (6.10) 
V -or 8 -a s J 

where A„ jS is a constant to be determined below. 

A n (z) is uniquely determined by these conditions, and can be written as 

A n (z) = £ A n , s ( a \ 1 ) (C(^ - 7.) + C(7.)) 



n-1 



, -IV 

+ E "-^dTa^Un-l-m + Un- (6.11) 



The coefficients A„ jS are determined by the linear equations 



E M a \ 1 I + E/n-i = (6.12) 



—al — a, 
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or, equivalently, 



^ ^n,s — —(U n -l)l2, 
s=l,2,3 

^ OZ s ^n,s = _ (fn-l)ll — (^n-l)22, 
s=l,2,3 

E «X- = -(^n-l)21. (6.13) 



s=l,2,3 



Of course, these linear equations are uniquely solvable as far as the genericity condition 
A 7^ is satisfied. 

Lastly, the hierarchy is defined by the system of Lax equations 

[d tn -A n (z), U(z)}=0 (6.14) 

for U (z) and the differential equations 

d tn t a 8 + t a a ^ 1 ) = « n , i t a i , (6.15) 

for t ot s = (a s , 1). The differential equations for j s reduce to 

^ 7 5 = -TrA n , s ( a \ 1 ] = (6.16) 

as expected. One can derive the zero-curvature equations [dt m — A m (z), dt n — A n (z)] = 
by the same procedure as in the previous case. 

6.3 Riemann-Hilbert problem and Grassmann variety 

The Riemann-Hilbert problem and the mapping to an infinite dimensional Grassmann 
variety can be derived in almost the same form as the previous case. The present case 
is conceptually rather simpler, because the poles 7 S do not move. To avoid confusion, 
we again move to the convention that the t-dependence is always explicitly indicated as 
A(t,z), A n (t,z), a s (t), etc. Note that 7 s 's are constant throughout the present setting. 

The Riemann-Hilbert pair consists of a Laurent series solution ip(t, z) and a global 
solution x(t, z) of the same auxiliary linear system. The former takes the form 

^(f, z) = <j>(t, z) exp(f] t n Jz~ n ), 0(t, z) = I + f] <f) n (t)z n . 

n=l n=l 

The prefactor <p(t, z) is connected with the generating function U (t, z) (the t-dependence 
is now shown explicitly) as U(t,z) = (f>(t, z)J<f>(t, z) -1 . The second solution x(t, z) of 
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the auxiliary linear system is characterized by the initial condition x(0, z) = I. One can 
prove, by the same technique as the previous case, that x(^ 5 z ) has essential singularity at 
z = and poles at z — 71, 72, 73, and behave as 



a s (0) 1 
-a s (0) 2 -a s (0) 



X (t,z) = \ x , s \ s \> ^ \( z - ls )- 1 + 0(l) (6.17) 



as z — > 7 S . These two solutions 2), ^(i, z) of the auxiliary linear system obeys a 
relation of the same form as (j4.13|) . On the other hand, since ^4 n (z)'s are trace-free, both 
x{t, z) and <f>(t, z) are now unimodular, i.e., 

det x(t, z) = det <j>(t, z) = 1. (6.18) 

Consequently, unlike the previous case, x(t, 2) has no degeneration point. 

We use the same Grassmann variety Gr to embed the hierarchy. The definition of the 
base point Wo (7, oc) for the present case, too, is essentially the same, except that we now 
use the three pairs (7 s ,a s ), s = 1,2,3, as the input. The basis {w n ,%j{ z ) \ n > 0, i,j = 
1, 2} of W (j, a) consists of the matrices w n ^j{z) of meromorphic functions on Y uniquely 
determined by these parameters as in the statement of Lemma E2 the third condition 
therein has to be modified as 

«•„,,( --.) A,,,,,. a \ 1 ) (z - ^y 1 + 0(1) [z^ ls ). (6.19) 
y -a z s -a s J 

By the correspondence 4>(t, z) 1— >■ W{t) = Wo (7, a(t))<f>(t, z), the hierarchy is converted to 
a dynamical system on the set M. of dressed vacua. The motion of W(t) again turns out 
to obey the same exponential law as (j5.9jl . 



7 Conclusion 

We have elucidated the status of the two elliptic analogues of the nonlinear Schrodinger 
hierarchy in the Grassmannian perspective of Sato and Segal and Wilson jT^j. Each 
of these systems are mapped to a dynamical system in the Grassmann variety Gr. The 
phase space of the dynamical system is the set M. of dressed vacua W = Wo (7, a)<p(z). 
The motion of the dressed vacuum W(t) = Wo (7(2), a(t))<p(t, z) under time evolutions of 
the hierarchy obeys a simple exponential law. This is just the restriction of universal expo- 
nential flows on the Grassmann variety itself. Thus the situation is fully parallel to many 
classical soliton equations that have been understood in the Grassmannian perspective. 
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It is straightforward to generalize the 2x2 system of the Krichever type to an r x r 
system [H]. In that case, one has to use several [/-matrices rather than a single one. The 
Tyurin parameters consist of r pairs (7^, ct s ) £ T x P r_1 , s = 1, . . . , r, of a point 7 S of T 
and an r dimensional directional vector ct s . As a special case, one can obtain an elliptic 
analogue of the so called N wave system, etc. 

If one does not insist on an explicit description of the system, one can generalize the 
results of this paper to an algebraic curve T of genus g with a marked point P . The Tyurin 
parameters for the construction of the Krichever type consist of 2g pairs (7^, a s ) £ T x P 1 
of a point of T and a two dimensional directional vector [Hj. The construction a la 
Enriquez and Rubtsov requires 3g, rather than 2g, pairs of Tyurin parameters |2j. Upon 
choosing a local parameter z in a neighborhood of Pq, one can start the construction of the 
fundamental matrix A(P) (P £ T) of meromorphic functions on T and the matrix U(z) 
of Laurent series. A convenient choice of z is to define it as the (multivalued) primitive 
function z(P) = Jp Q 10 of a holormophic differential u on T without zero at P . The 
matrices A n (P), like A(P), are characterized by a set of conditions on the poles. Namely, 
they are matrices of meromorphic functions on F with poles at Pq and 7 s 's, and behave 
as 



theorem. Formulating these systems in a more explicit form is a problem left for future 
research. 

Lastly, let us mention some other approaches to soliton equations associated with 
algebraic curves. Ben-Zvi and Frenkel pQ and Levin, Olshanetsky and Zotov ^U] propose 
to construct those equations as a 1 + 1 dimensional analogue of the Hitchin systems 
[3]. The framework of Ben-Zvi and Frenkel is conceptually similar to ours, though they 
use a Grassmann variety in a different way. The work of Li and Mulase |12t |§| is also 
closely related to the present issue. Our construction of dressed vacua has obviously 
a counterpart in their description of commutative rings of differential operators in the 
language of infinite dimensional Grassmann variety. 

Acknowledgements 

I would like to thank Takeshi Ikeda and Takashi Takebe for discussion and many com- 
ments. This work was partly supported by the Grant-in- Aid for Scientific Research (No. 



MP) = U(z(P))z(P)- n + 0(z(P)) (P -> Pq), 




The existence and the uniqueness of these matrices are ensured by the Riemann-Roch 



35 



14540172) from the Ministry of Education, Culture, Sports and Technology. 



References 

[1] D. Ben-Zvi and E. Frenkel, Spectral curves, opers and integrable systems, Publ. 
Math. Inst. Hautes Etudes Sci. 94 (2001), 87-159. 

[2] B. Enriquez and V. Rubtsov, Hecke-Tyurin parametrization of the Hitchin and KZB 



systems, e-print arXiv:math. AG/9911087 



[3] E. Frenkel, Five lectures on soliton equations, Surveys in Differential Geometry, vol. 
4, pp. 27-60 (International Press, 1998). 

[4] N. Hitchin, Stable bundles and integrable systems, Duke Math. J. 54 (1990), 91-114. 

[5] I.M. Krichever, Commutative rings of ordinary linear differential operators, Funct. 
Anal. Appl. 12 (1978), no. 3, 175-185. 

[6] I.M. Krichever and S.P. Novikov, Holomorphic vector bundles over Riemann surfaces 
and the Kadomtsev-Petviashvili equation. I, Funct. Anal. Appl. 12 (1978), no. 4, 
276-286. 

[7] I.M. Krichever and S.P. Novikov, Holomorphic bundles over algebraic curves, and 
nonlinear equations, Russian Math. Surveys 35 (1980), no. 6, 53-80. 

[8] I.M. Krichever, Vector bundles and Lax equations on algebraic curves, Commun. 
Math. Phys. 229 (2002), 229-269. 

[9] Y. Li and M. Mulase, Prym varieties and integrable systems, Commun. Anal. Geom. 
5 (1997), 279-332. 

[10] A.M. Levin, M.A. Olshanetsky and A. Zotov, Hitchin systems - symplectic Hecke 
correspondence and two-dimensional version, Commun. Math. Phys. 236 (2003), 93- 
133. 

[11] L.J. Mason and N.M.J. Woodhouse, Integrability, self-duality and twistor theory, 
Clarendon Press, Oxford, 1996. 

[12] M. Mulase, Category of vector bundles on algebraic curves and infinite dimensional 
Grassmannians, Intern. J. Math. 1 (1990), 293-342. 



36 



[13] E. Previato and G. Wilson, Vector bundles over curves and solutions of the KP equa- 
tions, Proc. Symp. Pure Math. vol. 49, part I, pp. 553-569 (American Mathematical 
Society, 1989). 

[14] M. Sato and Y. Sato, Soliton equations as dynamical systems on an infinite di- 
mensional Grassmannian manifold, Lecture Notes in Num. Appl. Anal., vol. 5, pp. 
259-271 (Kinokuniya, Tokyo, 1982). 

[15] G.B. Segal and G. Wilson, Loop groups and equations of KdV type, Publ. Math. 
IHES 61 (1985), 5-65. 

[16] A. Tyurin, Classification of vector bundles over an algebraic curve of arbitrary genus, 
AMS Translations II, Ser. 63, pp. 245-279. (American Mathematical Society, 1967). 

[17] V.E. Zakharov and A.V. Mikhailov, Method of the inverse scattering problem with 
spectral parameter on an algebraic curve, Func. Anal. Appl. 17 (1982), No. 4, 247- 
251. 



37 



